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Abstract 

We apply a Coupled Markov Chain approach to model rating transitions and thereby 
default probabilities of companies. We estimate parameters by applying a maximum likeli- 
hood estimation using a large set of historical ratings. Given the parameters the model can 
be used to simulate scenarios for joint rating changes of a set of companies, enabling the use 
of contemporary risk management techniques. We obtain scenarios for the payment streams 
generated by CDX contracts and portfolios of such contracts. This allows for assessing the 
risk of the current position held and design portfolios which are optimal relative to the risk 
preferences of the investor. 

1 Introduction 

The analysis and valuation of structured credit products like Collateralized Debt Obligations 
(CDO) and Credit Default Swap Indices (CDX) gained significant importance during the 
sub-prime mortgage crisis since 2007. Financial institutions still hold many products for 
which the risk exposure is unknown. 

The core problem of the recent sub-prime mortgage crisis and the connected turmoil on 
financial markets is that the risk of portfolios of certain asset backed securities were grossly 
underestimated. In particular the implicit underestimation of dependencies of default events 
and the resulting overestimation of the credit quality of products resulting from (repeated) 
tranching proved to be highly problematic. 

It turned out that common risk factors related to the U.S. economy in general and to 
the housing market in particular lead to huge losses for a certain type of financial products 
(mainly mortgage backed CDOs), that were previously rated super senior. 

To avoid such situations in the future it is crucial to model correlations between de- 
fault events and the common risk factors carefully and fully understand the effects of these 
dependencies on complicated structured products. 

There exist many models for describing the joint default behavior of companies in the 
literature (for an excellent survey and classification see for example [5] [6] [7] El [9] or [14) ) . In 
this paper we use a so called coupled Markov chain model (originally published in [11]) to 
capture dependency between credit events. Rating transitions of companies are modeled as 
discrete Markov chains and dependence between rating (and therefore also default) events 
is achieved by coupling the involved Markov chains. 
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Although rating based models have been criticized (see [Tl [4l Il3l 117] ) for their sluggish 
response to fast evolving events (like a rapid credit quality deterioration of a company), they 
have the advantage that they are based on a (usually) reasonable assessment of the credit 
worthiness of the debtor (combining a range of qualitative and quantitative approaches) and 
are therefore a signal for the credit quality of the debtor, which is consistent over time as well 
as among different debtors. On the other hand structural models like the KMV model (see 
for example |4]) that strongly depend on the asset value of the company under consideration 
suffer from the problem of overreaction to changes in the market - especially in time of crises 
or boom. This leads to gross over- or underestimation of credit risk associated with a loan 
portfolio. Rating based approaches offer a more stable assessment, since big rating agencies 
aim for ratings that reflect default probabilities through the cycle. 

Furthermore rating based evaluations are of increasing importance since pending new 
banking regulation uses ratings as an important input for calculating capital requirements 
for banks (see [3]). 

The most commonly used rating based approach to modeling credit risk is used in the 
popular CreditMetrics software. The main idea behind CreditMetrics is similar to the ideas 
in this paper: the current rating of a company influences the default probability in the 
next period. However, the crucial difference is that in the CreditMetrics model the ratings 
are treated merely as signals which are used to calibrate a multi-variate normal copula 
model for the joint rating transitions (see [14] for a discussion). Hence, the model assumes 
a parametric distributional model for the joint behavior of the rating changes. While this 
assumption might be reasonable in a normal economic environment, it often fails to model 
the default behavior appropriately in times of crisis (where the corresponding distributions 
are typically non-normal due to lack of symmetry and heavy-tailedness) , i.e. in situations 
when a correct assessment of downside risk is needed most. 

The coupled Markov model discussed in this paper avoids these assumptions and mod- 
els the joint behavior of the debtors in a non-parametric setting, constrained only by the 
assumed functional form of the rating transitions. The model is similar to the models pre- 
sented in |15l 1161 121] in the sense that it is a rating based model that can be represented as 
a mixture model. 

In |15l 116] the authors propose a generalized linear mixed model (GLMM) for rating 
transitions which is estimated using Bayesian techniques (GIBBS sampler). The model 
systematic risk factors as a combination of fixed and random effects and also allows for 
serial correlations in the unobserved risk factors and hence for so called rating momentum 
on a macroeconomic scale. 

In [21] a model for a continuous credit worthiness variable is proposed which are trans- 
lated to discrete ratings by identifying a rating class with an interval of the credit worthiness 
score. The continuous credit worthiness variable is allowed to depend on obligor specific as 
well as macroeconomic factors, the latter are used to model dependencies in rating transi- 
tions of different obligors. 

The contribution of the paper is twofold: 

1. We develop methods of statistically estimating the parameters of the credit migration 
model published in [11] and thereby enabling the use of the model for risk measurement 
and management purposes via Monte Carlo sampling from the fitted model. 

2. Returns from credit portfolios typically have highly skewed distributions - a character- 
istic further amplified by repeated tranching as observed in the markets in the recent 
years. One aim of this paper is to demonstrate that because of this fact risk manage- 
ment techniques are essential, when deciding on portfolios of structured credit products 
- not only ex-post to determine how much risk capital is needed to secure existing po- 
sitions but also ex-ante to design optimal portfolios relative to the risk preferences of 
the investor. 

In the second part of the paper we use Monte Carlo simulation techniques to obtain 
scenarios for the payouts of financial products like CDX contracts. These scenarios 
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are subsequently used to capture portfolio effects inherent in the considered structured 
loan products and their impact on the risk characteristics of portfolios containing 
such products. Using the distributional information obtained by the credit migration 
scenarios, we outline an optimization framework which can be used for active risk 
management for portfolios of CDX tranches and demonstrate the power of the approach 
by numerical examples. 

This paper is structured as follows: Section [2] is devoted to a short review of the coupled 
Markov chain approach. In Section [3] we discuss a maximum likelihood approach which 
is subsequently used to estimate the parameters of the model from empirical data. In 
Section [4] financial products based on Credit Default Indices are described while section [5] 
uses the coupled Markov chain model to calculate fair spreads for standard tranches of the 
iTraxx index. Finally in Section [6] portfolio optimization techniques are applied to portfolios 
of structured credit products - explicitly modeling the (joint) risk in the portfolios, while 
Section [7] concludes the paper. 

2 The Model 

This section gives an overview of the coupled Markov chain model introduced by [11] and is 
included to keep the paper self contained. The interested reader is referred to the original 
work for further details on the model. 

Consider a diversified portfolio consisting of debt obligations of different firms. The 
debtors are non-homogeneous in their credit ratings and they belong to different industry 
sectors. Assume that there are M non-default rating classes. The ratings are numbered 
in a descending order so that 1 corresponds to the safest class, while M is next to the 
default class. For example, in terms of the rating scheme of Standard and Poor's (S&P), 
1 < — ► AAA, 2 < — > AA, 3 < — ► A, 4 < — > BBB, 5 < — ► BB, 6 < — > B, 7 < — > CCC, 
8 < — > CC, 9 < — > C and 10 * — ► D with M = 9. 

The methodology suggested in [TT] describes joint rating transitions of companies with 
different rating classes belonging to different industry sectors, such that 

1. migrations of companies having the same credit rating are dependent; 

2. evolution of companies through credit ratings are dependent; 

3. every individual migration is governed by a Markovian matrix, which is the same for 
all the companies. 

The model assumes that the rating process of each debtor is Markov with identical rating 
transition matrix and that these processes are coupled in such a way that they are statisti- 
cally dependent through their common dependence on systematic factors. The specification 
of this dependency structure can be seen as the main contribution of the model. 

Note that there is some empirical evidence hinting to the fact that the Markov assumption 
of credit ratings does not always hold (see [T] 1131 I17j). Nevertheless this assumption is 
implicit in most of the popular credit risk models (also in the structural models) since it 
greatly simplifies the analysis. As the assumption does not seem to be too wrong (as shown 
in [T3]), we do not consider more complicated models. 

Let pij be the probability of transition within one year from the i-th credit rating to the 
j-ih. In particular, p it M+i is the probability that a debtor who is currently rated with i-th 
credit rating defaults by the end of the year. The M x (M + 1) transition matrix P = (pij) 
can be estimated using one of the various techniques proposed in literature for this purpose 
and is also reported by the rating companies themselves. Since we are mainly concerned 
with the coupling of the rating processes we assume the matrix P to be given. 

Consider X(t), t > 1, a discrete-time Markov chain evolving in the state space {1,2,..., M-\- 
1} governed by the probabilities p it j, whose absorbing state is M + 1. Suppose the evolution 
of the total portfolio is modeled by a multi-dimensional random process X 1 = (X{, . . . , X l N ) 
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whose components evolve like X(t), but are dependent. Denote by s(n) the industry sector 
of firm n. 

Since X l n is modeled as time-homogeneous Markov process, it would suffice to look at 
the transition from time t = f to time t = 2 and we could therefore drop the index t for the 
description of the probabilistic model. However, because we need the time index t in the 
portfolio selection models in Section [5] and [6] we will keep the subscript to have consistent 
notations throughout the paper. 

The starting point of the analysis for the proposed model is the commonly used division 
of risk factors in an idiosyncratic and a systematic factor^ 

To model dependencies between the rating transitions the discussed model also follows 
this intuition and splits the rating movement of company n at time t in two components: 
an idiosyncratic component r]„ and a systematic component £„. The distribution of both 
components depends on the current rating class of the company and the non-trivial joint 
distribution of the rf n is used to model dependency between the rating processes of the 
companies in the portfolio. 

The mixing between the two components is achieved via a Bernoulli variable S„ depending 
on the rating class of company n at time t — 1, i.e. as well as the industry sector s(n) 

of company n. We write 

P(<5^ = f | X*" 1 = m, s(n) = s) = q m , 3 
and organize the respective probabilities in the matrix 

Q = (?m,s)l<m<(M + l),l<s<S- 

Note that the sectoral classification used in this paper as an input for the model can be re- 
placed by any arbitrary (discrete) classification scheme for companies without structurally 
changing the model. Possible alternative classifications could for example be size or geo- 
graphic origin of the company. 

The rating of company n at time t is given bjj^] 

X = Sn£n + (1 - SnWn- 

The variables 8^ and are independent of everything else, while the variables rf n are 
independent of the 8^ and the {„ but have a non-trivial dependency structure within a given 
time period. The marginal distribution of the and the rj n is dependent on the current 
rating class of the company n and the respective distribution resulting from the transition 
matrix P, i.e. 

Hin = i)=P x l-\v Vie {f,...,M + l} 

and 

F(f) t n=i)=P xi -i ii , V«e{l,...,M + l}. (f) 

To describe the dependency structure of the if n we need several steps. First note that 
any transition of a discrete Markov chain rf n may be divided in two moves. One of them 
determines only the tendency of the move (i.e. up or down), while the other one determines 



1 Similar to the classical one factor Gaussian model where the asset value used to determine the probability of 
default is modeled in dependence of an idiosyncratic factor and a systematic factor which is common to all the 
companies 

2 Note that in |lll the ratings were described by the slightly different model 

= + (1 - 

where there are only (M + 1) variables n. Which variable is relevant for firm n was determined by the rating class 
of n in the previous time period. In this way dependency was introduced between the companies in the same 
rating class. Although the approach followed here uses a separate variable rj^ for every company, dependency 
between these variables in modeled in an analogous way (see discussion below). 
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the magnitude of the change. To be more specific let \m f° r fn = 1, . . . , M be Bernoulli 
tendency variables for the Markov process rf n . These variables determine whether a non- 
deteriorating move takes place for company n with -X^ _1 = m, i.e. 



t 



1, % > »7n 



t+1 



Xm \ 0, rf n <rf+\ 

In other words, xln is the indicator of a non-deteriorating move at time t. If the probability 
of transitions are given by the matrix P as in |TJ the tendency variable has a Bernoulli 
distribution whose probability of a success reads 

i 

Pi :=£w,i=P(Xn = l|^" 1 =0- 

2=1 

Observe that there are not TV but only M tendency variables per time period t, meaning 
that each of these determines the tendency of all the variables n^ for the subset of companies 
belonging to a specific rating class. This renders the variables r\ n dependent via their 
common components x = (x*> • ■ • > Xm) and the non-trivial joint distribution of the same. 

Conditioned on the tendency, the second move determines the magnitude of the change 
and follows the distribution 



^(Vn =j\X n =i, Xx*- 1 =1 ) = 

and 

nni=i I x'- 1 = i, x * t-i =o) 



Pij/Pt, j<h 
0, j > i 



Pi.j/^-Pi), 3>i 
0, 3 < i. 

The above deconstruction of a rating move in a tendency and a magnitude is useful, since the 
tendency part can be used to model the dependencies between different debtors via the joint 
distribution of x* == (x'i ■ • ■ > Xm)- Given the tendencies the distribution of the magnitude is 
completely determined by the marginal distribution of the rating transitions, i.e. the matrix 
P. Note that the variables rf n \ X x t-i ar e independent of one another. 

As mentioned above the presented model can be viewed as a Bernoulli mixing model, i.e. 
a model where the default probability of company n is dependent on some common factors 
X and 

P(X* = (M + l) | X )=Pn(x) 
with the influence function p n : range{\) [0, 1]. In the above case we have 

. Px^r 1 ,m+i 



Pn(x)=P X t n -\M+l<lx t n - 1 ,s(n) + t 1 ) i 1 ~ ^X*" 1 ,.(n)) 



P xi-' 



Unlike in other Bernoulli mixture models the common factors x are n °t explicitly modeled 
as economic factors but are mere statistical phenomena. 

In the original paper pairwise correlations between the tendency variables were assumed 
to be known and the joint distribution of these variables was found by an optimization 
approach (see [11] for further details). The probabilities in Q were manually tuned for 
small models, to demonstrate the models flexibility and expressiveness. However to use the 
model in practice numerically tractable estimation routines are needed to obtain realistic 
parameter estimates for Q and P x . 

For a model instance with M non-default rating classes and S industry sectors and known 
transition matrix P, there are S x M unknown parameters in Q and 2 M — (M + 1) for the 
specification of the joint distribution x* — (xL • • • > Xm)- I n the next section we will describe 
how to estimate these parameters given data on joint rating transitions via a maximum 
likelihood approach. 
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3 Estimation of Parameters 



We consider the model 

Xi = Sl£ + (1 - SDvl, (2) 

described in Section [2] where 

F(C n = i)=nv x i-\ i )=Px t -\i 

is assumed to be known. Given a set of realizations of the rating process for TV firms 
X = (x 1 ,...,x T ) £ R (JV,T) of the N dimensional process X 1 for different rating and industry 
classes over a period of T time steps, we want to estimate the parameters Q — {qi,j)i,j and the 
joint probability mass function P x for \ = (Xi> • ■ • > Xm) by a maximum likelihood approach. 
Since we model rating transitions via a coupled Markov chain approach, X provides us with 
N(T — 1) realizations of the process Q (for different industry classes and rating classes), 
but only with (T— 1) (hidden) joint realizations of \ — (Xii • ■ • > Xm). It follows that we can 
analyze every time step separately, but have to take into consideration the joint behavior of 
moves within each step. In the following we derive the likelihood function by conditioning. 
The likelihood of the matrix X given the matrices Q and P x , is 

T 

L(X; Q, P x ) = P (X 1 = x* | X'- 1 = x'- 1 ) (3) 

i=2 

and by the law of total probability 

where X 1 — (X{, . . . , X l N ) and x* = (Xii ■ ■ ■ i x\i)- We further have because of independence 
of x* and rf 

F(X* = a?*|x* = X,^"*" 1 = x' 1 ) = 

FiX* = x'lx* =X,S* = 5, X*- 1 = x 4 - 1 )?^* = 6) 

Je{o,i} JV 

with 5* = (#,...,<&). 

To calculate the above sum it is necessary to divide the companies into groups and 
consider all possible realizations in the considered subgroups separately. Define I 1 (s, mi , 7712) 
to be the number of companies in sector s which move from class mi to class mi in period 
t. We start by analyzing these subgroups. 

First fix a t and (s, mi, mi) with mi < mi and Xmx — 1- We calculate the probability 
that the I*(s, mi, mi) companies can move from mi to mi as a function of the parameters to 
be estimated. Since the tendency variable x mi = 1, the only possibility for a deterioration 
mi — * mi to happen is that <5„ = 1 for all the corresponding companies (all other realizations 
of 8^ yield probability 0). The probability is therefore 

I t (s,m 1 ,m 2 ) I t (a,m 1 ,m 2 ) 
™i,s Pmi,m2 

Now consider the case that mi > mi and X mi = 1. In this case the corresponding 
companies have the possibility to move from mi — * mi either by a realization of f„ (if 
the Sn — 1) or by a realization of r)„ (if 8 n =0). Since all the combinations of these two 
possibilities have to be considered, the probability is 

ft jt 

/ I \ ( (1 - q mi ,s)p mi ,,n 2 \ n , ,Z*- n _ I* / QmusiPmi ~ I) + 1 \ 

/ j I _ I I + I Wmi, sPm 1 ,m 2 ) — Pm 1 ,m 2 \ i ] 

^ V 71 / V Pmj / V Pm! / 



G 



where we abbreviate P (s , mi , m 2 ) by I . 

We now split up all the moves in a period t according to industry sector, starting class 
and rating class in (t + 1) and by the above argument get the following formula 



P(X t = x i \ X t =x,X t - 1 =x*- 1 



S M M+l 

=n n n /<y~ i > s > m i> m 2>; ( 2'- p x 

s = l mi = l T712 — 1 



with 



f(x l ,s,m 1 ,m 2 ,;Q,P x 



Pm 1 , 



Pin 1 , 



Pm, 



, 9mi ,sPmi ,7712 ' 



mi > m 2 , Xmi = 1 

mi < m 2 , Xmi = 
otherwise. 



The likelihood function (|3| can now be written as 



L(X;Q,P X ) = II E Px(x t = x)nX* = x t \ X t = X,X i - l = x t - 1 ) 
t=2 xe{o,i} M 

T 

= II E P x(x' = X) II f(x t ~ 1 ,s,m l ,m 2 ,;Q,P x 

t — 2^^10,1}^ s,mi,m 2 



The above function is clearly non-convex and since it consists of mix of sums and products 
this problem also cannot be solved by a logarithmic transform. 

Maximizing the above likelihood for given data X in the parameters P x and Q amounts 
to solving the following constrained optimization problem 



max Q ,p x L'(X;Q,P X ) 

s.t. q m:S £ [0, 1] 

Ex:x«=i P x(x) =P^i> Vi : 1 < i < M 

Ex:xi=o P x(x) =1-Pt- 



where 



(4) 



L'(X;Q,P X ) =^log\ ^ P x (x* = x) Yl 1 , s,m 1 ,m 2 , ;Q, P x 

t — 2 \-^£{o,l} M s,m 1 ,m 2 



and 



' / « mi ,s(pt 1 - 1 )+ 1 \ lt 



f(x ,s,m 1 ,m 2 ,;Q,P x ) = < 



gm 1|3 (P mi -l) + l 

Pm, 



mi > m 2 , Xmi = 1 

mi < m 2 , Xmi = 
otherwise. 



The modified log-likelihood L' is numerically more tractable than the original likelihood 
function, because of the logarithm and the omission of the terms p mi>m2 . Furthermore it is 
obvious that both extremes are at the same value and therefore yield the same maximum 
likelihood estimator. 

Since the problem is non-convex and the amount of parameters is too large to employ 
standard global optimization techniques, we use heuristic global optimization techniques to 
solve Q. A detailed description of the algorithms used and their properties can be found 
in \m. 
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For the rest of the paper we will assume that the parameters (Q,P X ) axe estimated via 
a maximum likelihood approach. The estimation is based on yearly historical rating data 
for 10166 companies from all over the world produced by Standard & Poor's over a time 
horizon of 23 years (1985 - 2007). The data set comprises of 87296 observations of rating 
transitions (not all the companies are rated over the whole time horizon). We number the 
10 Standard & Poor's rating classes from 1 to 10, 1 being AAA and 10 the default class. To 
reduce the number of parameters the 10 rating classes are clubbed in the following way: 1 
<- AAA, AA, 2 <- A, 3 <- BBB, 4 <- BB, B, 5 <- CCC, CC, C, 6 <- D to obtain a reduced 
set of 6 rating classes (M — 5). We estimate the matrix of transition probabilities P by 
simple counting as 



P = 



0.9191 
0.0335 
0.0080 
0.0039 
0.0023 




0.0753 
0.8958 
0.0674 
0.0092 
0.0034 




0.0044 
0.0657 
0.8554 
0.0794 
0.0045 




0.0009 
0.0036 
0.0665 
0.8678 
0.1759 




0.0001 
0.0006 
0.0011 
0.0244 
0.6009 




0.0001 
0.0009 
0.0016 
0.0153 
0.2131 
1 



Note that in our data set we have transitions from every (non default) rating class to every 
other rating class and therefore have no problems with probabilities which are estimated to 
be equal to 0. 

Additionally we incorporate information on industry sectors by considering the following 
6 industry classes (according to SIC classification scheme) Mining and Construction (1), 
Manufacturing (2), Transportation, Technology, and Utility (3), Trade (4), Finance (5), and 
Services (6). 

The choices above leave us with 62 variables to be estimated. The estimated matrix Q 
is given by: 



0.1981 
0.3854.10" 7 
0.2732. 10" 10 
0.0299 
0.0816 



0.0818 

0.2008.10" 

0.3337.10" 

0.0487 

0.1739 



0.0138 
0.0655 
0.0344 
0.1518 
0.4437 



0.0001 

0.8457.10" 

0.0005 

0.0341 

0.4092 



0.1467 
0.0609 
0.0494 
0.0076 
0.2482.10" 



0.3089 

0.0356 

0.1752.10" 

0.0300 

0.0000 



while the joint probability function for the x can De found in Table [T] 



4 Credit Default Swap Indices 

Below we describe how to use the econometric model for credit rating transitions to design 
portfolios of structured credit products taking the risk resulting from the uncertainty of fu- 
ture payoffs into account. More specifically we want to construct portfolios of Collateralized 
Debt Obligations (CDO) based on Credit Default Swap Indices (CDX). 

In this section we outline the setting for the portfolio optimization carried out in Section 
|6j discuss the concept of a CDX, and specify the payoff function of a CDO based on a CDX, 
which in turn is based on the credit events of single companies. 

The payoffs of CDX contracts are based on the performance of portfolios of so called 
Credit Defaults Swaps (CDS) contracts. 

The basic idea of a CDS contract is that the buyer of the contract transfers her credit 
risk associated with the loans to a certain third party (usually called the reference entity) to 
the seller of the contract, i.e. the seller insures the buyer against the possible losses resulting 
from a credit event of the reference entity. In case of a credit event the losses of the buyer 
are covered by the seller of the contract (cash settlement^] A credit event can be a default 

3 There are also other forms of settlement for CDS contracts other than the cash settlement. Some contracts 
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Rating Class 



Rating Class 



1 2 3 4 5 



Probability 



1 2 3 4 5 



Probability 





1 

1 

1 1 

10 

10 1 

110 

111 

10 

10 1 

10 10 

10 11 

110 

110 1 

1110 

1111 



0.0058 
0.0001 
0.0001 



0.2670e-004 

0.2743e-005 

0.4057e-006 

0.9493e-004 

0.0335 

0.3368e-004 

0.1811e-004 

0.0001 

0.2545e-006 



0.0081 
0.0005 
0.0312 
0.1335 



1 
1 1 
10 10 
10 11 
10 10 
10 10 1 
10 110 
10 111 
110 
110 1 
110 10 
110 11 
1110 
1110 1 
11110 

11111 



0.4814c-005 
0.2059e-004 
0.1526e-004 
0.2533e-005 
0.2776e-005 
0.9985e-006 
0.4697e-004 
0.1352e-005 
0.4803e-004 
0.0219 



0.0002 
0.0409 
0.0341 
0.0001 
0.0011 
0.6885 



Table 1: Probability function of x, i-e. joint distributions of up and down moves in the different 
rating classes. The first 5 columns of the above table describe the event, i.e. the combination 
of up and down moves in the different classes (1 being a non-deteriorating move, while is a 
deteriorating move), while the last column gives the estimated probability. 

of the reference entity or any other event decreasing the credit quality such as failure to pay 
or restructuring. We will focus on default events as the only credit events in this paper. 

In exchange for the insurance provided by the seller of the contract the buyer pays a 
certain percentage of the insured amount (the notional) every quarter. Note that a CDS 
contract often takes the form of a mere bet: the two parties engage in the contract without 
the existence of underlying loans to the reference identity. 

A CDX is a portfolio of simple CDS contracts. In all cases observed in practice the 
weighting of the individual contracts is uniform, i.e. every CDS contract contributes equally 
to the overall notional of the CDX contract. As of now there are two companies providing 
CDX contracts (differing in composition, maturity as well as date of issue): iTraxx for 
European companies and DJ CDX for the U.S. market. 

In this paper we analyze so called CDX tranches, which are synthetic collateralized debt 
obligations based on CDX. Each buyer of a CDX tranche buys a particular segment of the 
CDX loss distribution. Typical tranches are equity (0-3%), mezzanine (3-6%), senior (6-9% 
and 9-12%) and super senior (12-21%). If the tranching is designed like above the first 
three percent of losses are absorbed by the equity tranche, the next three percent by the 
mezzanine tranche and so on, i.e. if for example 2% of the notational value of the CDX 
defaults, then only the buyers of the equity tranche face losses. The buyers of the different 
tranches receive different spreads on their invested capital. 

We assume that we have an asset universe of N = {1, . . . , N} C N companies which can 
be used as reference entities of CDS contracts included in the CDX. The payment streams 
generated by the contract will be viewed from the risk buyers side, i.e. positive numbers 
indicate profits for the risk buyer of the CDX tranche. 

specify that in case of a credit event of the reference identity the buyer of the CDS hands the seller a bond with 
face value equal to the notional amount of the CDS contract in exchange for a cash payment of the notional 
(physical settlement). For simplicity we will stick to cash settlement in this paper. 
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To describe a specific tranche C of a CDX contract we have to know the following 
parameters: 

1. The assets Ac C J\f that make up the CDX and the notional value N(C,t,n) with 
which the CDS corresponding to asset n £ Ac enters the CDX. Note that the notional 
value of a company n can drop to zero due to a credit event during the period for 
which the CDX is issued, hence the time dependence of the notional. We denote the 
notional value of the whole CDX at time t by N(C,t). 

2. The tranching specified by the so called attachment point ac G [0, 1) and the detach- 
ment point /3c £ (0,1], i.e. the point where the contract begins absorbing losses and 
the point where everything in the tranche is lost. 

3. The maturity Tc of the contract. 

4. The spread Sc of the contract, i.e. the positive difference of paid interest rates to the 
risk free interest. 

5. The recovery rate Rc, i.e. the fraction of the notional assigned to a CDS which can 
be recovered in the case of the default of the corresponding reference entity. For the 
sake of simplicity we assume this quantity to be non-stochastic and the same for all 
companies n and note that an extension to non-uniform recovery rates would require 
only trivial modifications of the methods presented in this paper. Even the most 
realistic case - namely stochastic recovery rates - could easily be incorporated into the 
proposed Monte-Carlo framework (see for example |19j for a non-parametric model for 
default rates). 

Hence we can specify our contract as C = (Ac, &c, 0c, Sc,Tc, Rc)- 

To simplify the problem, we assume that the contributions of every individual debtor to 
the CDX are the same, i.e. 

As mentioned earlier, in practice this is not a restriction, since all CDX products currently 
available on the market have this simple structure. It will become obvious that the equal 
weights are not crucial in our approach and all the evaluation formulas can be extended in 
a straightforward manner to the case where the assets have different weights. Having made 
this assumption the percentage of overall losses of the CDX contract that C belongs to can 
be written as 

J-in ~ 

t 

period t. We now can write N(C, i) as 



\M 

where D c is the number of defaults of reference entities of CDS contracts in the CDX until 



N(C,t) = < 



N(C,0)(f3 c 



a c < L* c < I3 C 



Pc— a c 

N(C,0), L* c <aa 
0, L* > Pc- 



As mentioned above it is furthermore assumed, that the recovery rate Rc is non- 
stochastic and the same for all the debtors in the CDX contract. This assumption simplifies 
the valuation, since the order of the default events does not influence the amount to be paid 
by a particular risk buyer. 

If we divide the time period from the starting of the contract until Tc into T equal time 
periods (usually years), the random discounted returns arising from the CDX trance C can 
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be written as 



T 



TZ(C) 



N(C, 0) + J2 r r [N(C, t)(S c + it) + R(N(C, 0) - N(C, t))it] 



t=i 



+ r T 



\N(C, T) + R(N(C, T) - N(C, 0))] 



where i t is 1 plus the risk-free rate of return in period t and rt is the discount rate for period 
t. 71(C) is random since it depends on the number of default events D c of the companies 
the individual CDS in the CDX are based on. Note that if a company defaults the recovery 
amount is kept by the issuer of the CDX and paid out at the maturity of the contract - in 
the meantime the capital yields the risk free interest for the risk buyer. 

The above payout structure corresponds to the so called funded version of CDX products, 
meaning that the notional amount has to be actually paid in the beginning and the capital 
yields an interest rate of (Sc + it) in period t. Alternatively one could consider unfunded 
CDX products where no payment is made at the beginning of the contract, but only when 
defaults actually occur. If there is no spread in the market for risk free capital the two 
contracts yield the same payouts. Since we use a simulation based pricing approach the 
methods presented here are not specific to a particular form of contract. 

Other than the losses from defaults of companies, the CDX tranche can also change value 
because of changes of creditworthiness of debtors that do not involve a default event. Such 
changes in price are not important if the contract is held until maturitjj^] If however the 
contract is sold before maturity these changes also have to be considered. In this paper we 
will assume that the buyer of the CDX tranche holds the contract until maturity. 



Usually CDS and thereby CDX tranches are priced by comparing expected values of the 
payoffs and the obligations from the contract. Using "risk free" evaluation techniques the 
expected values (under the risk free martingale measure) of these two quantities are usually 
equated to price the contract. Obviously the approach taken here does not permit analytical 
solutions in dependence of model parameters. In this section - as a first test of the presented 
model - we compute rates for different CDX tranches by equating payoffs and obligations 
as described above. 

We try to calculate "fair spreads" for CDO's on CDX tranches using the real world (or 
physical) probability measure obtained by samples of the model discussed in the previous 
section and compare the obtained spreads with the spreads offered at the market. This 
serves as a first empirical test of the model. More specifically, we use the model from section 
[2] and the parameters found in section [3] to generate scenarios for the rating transitions of 
all the reference entities of CDS products which are included in the CDX of interest. As 
a sampling procedure we use simple Monte Carlo, i.e. we simulate the randomness in the 
model to obtain equally weighted scenarios. Given these scenarios we can calculate scenarios 
for the return of a CDX tranche C by evaluating 1Z(C,(?) for the specific scenarios £ l of 
joint rating transitions of the companies in Ac- 

To check our model for consistency with market data, we calculate spreads for different 
tranches C = (Ac, etc , Pc , Sc , Tb, Re) of a given CDX product and compare these spreads 
with spreads offered in the market. More specifically we compute Sc such that 



4 This statement is not quite correct if we also consider the effects of changed book values due to rating changes 
of the reference entities making up the CDX contracts. This may result in changed tax burden due to mark to 
market valuation of the portfolio. For the sake of simplicity we will not considers such effects, since these also 
depend on the overall situation of the company managing the portfolio. 



5 Calculation of Spreads 



E(K(C*)) = 



(5) 
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Table 2: iTraxx Europe Series 9 spreads and upfront payments on 31.03.2008 relative to 3 month 
EURIBOR of 4.6 %. 



0%-3% 

3%-6% 

6%-9% 

9%-12% 

12%-22% 



5 Years 

-0.0415 (0.0915) 
0.0159 (0.0321) 
0.0129 (0.0180) 
0.0126 (0.0089) 
0.0126 (-0.0017) 



7 Years 

-0.0238 (0.0738) 
0.0175 (0.0388) 
0.0107 (0.0245) 
0.0099 (0.0138) 
0.0099 (0.0021) 



10 Years 
-0.0113 (0.0613) 
0.0255 (0.0424) 
0.0105 (0.0292) 
0.0082 (0.0178) 
0.0078 (0.0056) 



Table 3: Spreads that fulfill ^ for the different tranches and maturities of iTraxx Europe Series 
9. Differences to observed market spreads from Table [2] are given in brackets. 



for a CDX tranche C* = (Ac, otc, Pc, S c ,Tc, Rc)- Given a set of return scenarios and the 
other parameters of the contract C the spread can easily be found by simple line search 
methods. 

For our experiment we consider the iTraxx Europe Series 9 Version 1 for the maturities 
5 years, 7 years and 10 years. The iTraxx Europe consists of 125 companies with ratings 
varying from AAA to BBB (i.e. investment grade) at the initiation of an index series. All 
companies enter the iTraxx with the same weight. The spreads observed on the market are 
given m Table [2] 

For the sake of simplicity we adopt the iTraxx standard estimate for the recovery rate 
for all the companies and set R — 0.4. We furthermore assume a discount rate of 5%. Using 
this information we simulate 10000 credit rating scenarios over a horizon of 10 years for 
the 125 involved companies. Table [3] gives the fair spreads which fulfill |5| for the different 
tranches and maturities. 

The spreads for the respective equity tranches are computed under the assumption of 
prepayments given in Table [2]- this also explains the possibility of negative interest rates for 
these tranches. It can be observed that the calculated spreads are lower than the observed 
market spreads. This is intuitively quite clear, since we use the physical probability measure 
obtained by sampling from the coupled Markov chain model and calculate mere expectations 
under this model. Thereby the risk inherent in the products is not taken into account at all, 
and the fair spreads are lower than the actual market spreads. This is also reflected in the 



12 




Figure 1: Loss distributions of the mezzanine (a), the two senior (b) & (c) and the super senior 
(d) tranche of iTraxx Europe, Series 9 with 10 years maturity. The average loss is depicted by a 
solid line the VaRo.g by a dashed line and the VaR 9g by a dotted line. For (a) and (b) a kernel 
density estimator was used to depict the distributions. All the computations are performed with 
data till the 31.03.2008. 

differences to the market rates, which are consistently the highest for the equity tranches 
(i.e. the most risky ones) and gradually decrease with the seniority of the tranche - i.e. 
the observed differences to the market are due to varying riskiness of the products and can 
therefore be viewed as risk premiums paid by the market. 

As a byproduct of the above analysis we obtain the loss distributions of all the CDX 
tranches under investigation. The loss distribution of the 10 Year mezzanine, the 10 year 
senior tranches as well as the 10 year super senior tranches of the iTraxx are plotted in figure 
[l] The tranches are chosen to demonstrate the diminishing riskiness of the tranches with 
higher seniority: while the mezzanine and the first senior tranche still exhibit a long right 
tail indicating a high risk of big losses the distributions of the second senior and the super 
senior tranches are fairly concentrated around the mean with only very few loss scenarios. 
As mentioned above, this fact serves as an intuitive justification of the differences of the risk 
free spreads calculated in this section to the actual market rates. 

In the next section we will adopt a more explicit notion of risk and create decision 
models which can be used to design portfolios of CDX tranches with more favorable risk 
characteristics as the single contracts investigated in this section. 

6 Portfolio Optimization 

In this section we want to make use of the distributional information provided in the model 
(especially the dependency structures between different companies and thereby different 
CDX products). For the assessment of additional expected profits incurred by a financial 
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product relative to an existing portfolio the joint distribution of the returns is obviously 
irrelevant, since the expectation is a linear operator. However, the risk associated with 
a particular portfolio of assets is dependent on the exposure of the whole portfolio. This 
exposure in turn hinges on the dependencies between the positions in the portfolio. 

In practice CDX tranches and similar contracts are often analyzed in isolation without 
explicitly considering the portfolio of contracts which the newly acquired contract is part 
of. Our aim in this section is to design a portfolio wide risk assessment and optimization 
framework by using the available information on the joint default probabilities and thereby 
on joint returns of different CDX tranches. This information is contained in the scenarios 
used in the previous section to calculate fair spreads. 

We will demonstrate numerical examples which show that the portfolio view is clearly 
advantageous in the sense that it allows for substantial reduction of risk of the positions 
(holding returns fixed). 

For our numerical example we use the Conditional Value-at-Risk at level a (CVaR a ) as 
a risk measure. Given a random variable X, the Conditional Value-at-Risk at level a of X 
is defined as (see [2] and [20] ) 



CVaR a (X) = — — / F x 1 {x)dx 
1 - " J a 



= max la — —'E([X — a] ) : a G 
I a 

where Fx is the distribution of X and F x (x) = inf {t : Fx (t) > x} is its generalized inverse. 
In principle there is a wide range of available risk measures which can be effectively used in 
optimization frameworks. Our choice is motivated by: 

1. The favorable theoretical properties of CVaR. Since CVaR is a convex risk measure 
it leads to sensible decisions from a an economic point of view, favoring diversification 
over concentration - a property that for example the Value-at-Risk lacks (see |18|). 

2. The fact that CVaR is piecewise linear (in the finite scenario setting) making problem 
|6j| tractable numerically (see formulation |7]l below). 

3. The fact that CVaR is closely related to VaR which plays an important role in the 
Basel accord - the current regulatory framework for banks. In fact it is easy to see 
that CVaR is always bigger than VaR, making portfolios which have favorable CVaR- 
characteristics also attractive from a VaR perspective. 

We assume that there exists a range of different CDX contracts C = (Ci, . . . , Cd). The 
aim is to solve the following mean risk problem using the CDX tranches as investment 
possibilities. 

min m CVaR a (w T 11(C)) 

s.t. E(w T TZ(C)) > n , 

J2 W i = 1 

Wi E [a, b], Vi : 1 < i < d 

where by slight abuse of notation w T 1Z(C) = X^=i U)ilZ(Ci) and the decision variables are 
the portfolio weights w — (wi, . . . , Wd) G K d constrained by the bounds a, b G R. 

We know that in a scenario framework the above problem can be computed by solving 
the following linear optimization problem (see [20]) 

mim„ a + jkt E-^ z ° 



|s| ^ses ■ 

s -t- ji|E seSc ™ T r >m 

E™; = 1 

Wi 6 [a, b], Vi : 1 < i < d 

w T (, s -a <z s , VsG5 

z s > 0, 



(7) 
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Figure 2: Return distribution of the mezzanine tranche of iTraxx Europe, Series 9 with 10 years 
maturity. All the computations are performed with data until the 31.03.2008. 
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CVaRo.g 


VaRo.99 


CVaRo.99 


10Y / Mezzanine 


0.1599 


0.0333 


0.4240 


0.3197 


Run 1 


-0.1825 


-0.3009 


0.2835 


0.0375 


Run 2 


-0.0925 


-0.1299 


-0.0059 


-0.0522 



Table 4: The table shows the Value-at-Risk as well as the Conditional Value-at-Risk at level 
a = 0.9 and a — 0.99 for the loss distribution of the original CDO tranche as well as for the 
portfolio resulting from optimization run 1 (minimize CVaRo.g) an d from optimization run 2 
(minimize CVaRg.gg). 



where S is the finite set of equally probable scenarios, while z 3 are auxiliary variables 
required for the modeling of CVaR as linear program and £ s is the vector of losses from the 
contracts Ci in scenario s, i.e. £ s = — (7?.(Ci, s), . . . , R(G\, sj). 

Our setup for the numerical tests is the following: we choose one iTraxx tranche (specifi- 
cally the 10 year mezzanine tranche) and investigate its return distribution with the support 
of 10.000 simulated scenarios from our model with the estimated parameters. Figure [2] shows 
the distribution of the returns of this tranche. The distribution exhibits the typical long 
left tail that is commonly observed from portfolios of loans and synthetic products based on 
loan portfolios. While the mean discounted return of the tranche is 0.3133, the VaRo.Q and 
the VaRo.gg are considerably larger indicating the high risk associated with the product. 

The idea is to try to find a portfolio of all CDO tranches based on the iTraxx Europe, 
Series 9, which has the same expected return as the 10 year mezzanine trance, but more 
favorable risk characteristics (in terms of CVaR). We therefore solve problem |7| with the 
parameters /i = 0.3133, a* = bi = 0.5 for all 1 < i < d and a = 0.1 and obtain the portfolio 
shown in Figure [3] 

The optimal portfolio is clearly less risky (yielding the same expected return): while the 
distribution still is heavy tailed (on the left), it contains much less probability mass which 
is also reflected by the drastic improvements of CVaRo.g as well as VaRo.g relative to the 
original return distribution (see Table [4] for an overview). 

The optimization run was repeated for the parameter a — 0.99. The portfolio and loss 
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Figure 3: Part (a) of the figure shows the optimal portfolio found by solving problem ^ with 
parameters /i = 0.3133 and a = 0.9. The y-axis labels give the product (1/5Y is the equity 
tranche with 5 year maturity) while the x-axis depicts the portfolio weights. A bar to the left 
means that the portfolio is short in the respective position while a bar to the right indicates a 
long position. 



distribution for this run are depicted in Figure [4] The lower value of a is reflected in a 
structurally different shape of the return distribution. The portfolio yields a surprisingly 
low VaRo.g and VaRo.gg while still yielding the same expected discounted returns as the 
original CDX tranche. 

Since the risk characteristics of both portfolios are clearly preferable to the original one, 
it is shown that a portfolio view which explicitly incorporates dependencies between the 
products and the risk originating from the joint distributions of the losses yields better 
decisions and avoids positions that have high concentration risk. 



7 Conclusion 

In this paper we presented methods to estimate the parameters of the rating transition model 
proposed in [11) using a maximum likelihood approach. Subsequently the estimates are used 
to generate a set of scenarios for the joint rating transitions of the companies included in 
the iTraxx, Series 9 index. 

The calibrated model was used to show that even financial products, which are buffered 
against losses by tranching, can have quite unfavorable return distributions. This result 
hinges on the explicit modeling of the joint distribution of defaults, which is such that it is 
capable of capturing effects which resemble a crisis situation: i.e. systematic simultaneous 
defaults of many debtors. 

The developments on financial markets in the two years after the emergence of the sub- 
prime crisis in 2007 demonstrated that structured products like CDX tranches should not 
be viewed as stand alone instruments, which are constructed in a way that further risk risk 
control is superfluous but merely as single assets whose analysis has to be integrated in a 
company wide risk management. We used a bottom up approach to assess the riskiness of 
such structured products and thereby enabling a correct assessment of the risk inherent in 
these contracts. The presented framework can be used to derive strategies to mitigate the 
risk emanating from CDX tranches by solving a portfolio selection problem where the asset 
universe consists of various CDX contracts. It turned out that even in this simple framework 
with few hedging instruments a carefully selected portfolio of contracts clearly outperforms 
a single contract. 

The presented optimization methods offer a convenient framework which has two advan- 
tages: firstly it allows the use of arbitrary scenarios for company default as an input, which 
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Figure 4: Part (a) of the figure shows the optimal portfolio found by solving problem ([6]) with 
parameters fi = 0.3133 and a — 0.99. The y-axis labels give the product (1/5Y is the equity 
tranche with 5 year maturity) while the x-axis depicts the portfolio weights. A bar to the left 
means that the portfolio is short in the respective position while a bar to the right indicates a 
long position. Figure (b) shows the loss distribution of the optimal portfolio as well as the mean 
(solid line) the VaRo.g (dashed line) and the VaRo.gg (dotted line). 



makes it possible to use sophisticated statistical models to predict the joint rating (and 
default) behavior of companies and secondly it is flexible enough to reproduce the complex 
contract structures of CDX products without any problems. The flexibility of the model 
makes it possible to relax most of the assumptions made in this paper to make the model 
even more realistic. In particular it would be easily possible to incorporate stochastic mod- 
els for the interest rates and the loss given default as well as transaction costs and varying 
contract details. 

The mentioned flexibility of the approach together with the numerical feasibility of large 
problem instances make the outlined methods interesting for companies and investors who 
are faced with the problem of risk management in the context of corporate bonds and 
structures credit products. 
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